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Abstract 

In  this  paper  a  nonlinear  stress-based  controller  is  designed 
to  position  a  single-degree-of-freedom  shape  memory  alloy 
(SMA)  actuated  manipulator.  A  three  part  model  was  con¬ 
structed  based  on  the  dynamics/kinematics  of  the  arm,  the 
thermomechanical  behavior  of  SMA’s,  and  a  heat  transfer 
model  consisting  of  electrical  heating  and  natural  convection. 
A  variable  stracture  controller  is  designed  to  calculate  a 
desired  stress,  based  on  position  error.  The  desired  stress  is 
compared  with  the  actual  stress  which  is  computed  using  a 
Kalman  filter.  The  stress  error  is  then  used  for  control  via  a 
proportional-integral  (PI)  controller.  Numerical  simulations  are 
performed  to  investigate  tracking  performance  and  robustness. 
The  results  demonstrate  that  the  controller  design  is  highly 
accurate  in  tracking  both  stationary  and  time-varying  input 
signals. 

I.  Introduction 

The  complex  nonlinear  behavior  of  shape  memory  alloys- 
specifically  the  co-dependent  relationship  between  stress, 
martensite  fraction,  and  transformation  temperature-presents 
many  challenges  with  regards  to  control  design  for  SMA 
actuated  manipulators.  Some  methods  that  have  been  inves¬ 
tigated  include  linear  control,  pulse  width  modulation,  and 
nonlinear  control  [1]-[13].  The  shared  feature  in  each  of  these 
control  designs  is  that  position  is  used  as  the  feedback  state. 
In  this  work  we  present  an  alternative  solution  to  the  control 
problems  of  SMA  actuators  by  tracking  sttess  as  opposed  to 
position.  The  benefit  of  tracking  stress  in  rotating  systems  is 
that  the  nonlinearities  due  to  angular  position  can,  in  effect, 
be  cancelled. 

Feedback  linearization  with  the  addition  of  a  sliding  surface 
(i.e.,  a  variable  structure  controller)  is  used  to  calculate  a 
desired  stress  in  the  wire,  based  on  the  position  error,  and 
is  compared  with  the  actual  stress  that  is  computed  using  a 
Kalman  filter.  The  stress  eiTor  is  then  used  for  feedback  control 
via  a  proportional-integral  (PI)  controller. 

To  investigate  stress-based  control  performance,  a  three  part 
model  was  constructed  in  MATLAB/Simulink  based  on  the 
dynamics  of  single-degree-of  freedom  SMA-actuated  arm.  The 
model  incorporates  the  dynamics/kinematics  of  the  arm,  the 
thermomechanical  behavior  of  SMA’s,  and  a  heat  ttansfer 
model  consisting  of  electrical  heating  and  natural  convection. 
Based  on  the  model,  numerical  simulations  are  performed  to 
investigate  tracking  performance  and  robustness. 


Fig.  1.  Diagram  of  single-degree-of-freedom  SMA-actuated  arm 

II.  SMA  Manipulator  Model 
The  SMA  actuated  manipulator  model  consists  of  three 
parts:  (l)the  kinematic/dynamic  model,  (2)  the  SMA  ther¬ 
momechanical  model,  and  (3)  the  heat  transfer  model.  As 
depicted  in  Figure  1,  the  arm  is  actuated  by  a  bias  type 
actuator  constructed  with  SMA  wire,  pulleys,  and  a  linear 
spring.  Heating  the  SMA  wire  induces  a  negative  sttain  in  the 
material  creating  a  positive  torque  and  the  arm  rotates  coun¬ 
terclockwise.  Conversely,  cooling  the  wire  causes  the  arm  to 
rotate  in  a  clockwise  direction.  The  following  sections  present 
the  dynamic  model  used  in  subsequent  numerical  simulations. 
The  accuracy  of  this  model  has  been  demonstrated  in  the 
experimental  work  performed  by  Elahinia  and  Ashrafiuon  [9]. 

A.  Kinematic/Dynamic  Model 

The  dynamic  model  of  the  arm  including  spring  and  payload 
effects  is  represented  by 

Ie9  +  c6+[Tg{9)+Ts{6)]=Tyj{a)  (1) 

where  rw,Tg,  and  Tg  are  the  resulting  torques  from  the  SMA 
wire,  gravitational  loads,  and  the  bias  spring,  respectively,  and 
a  is  the  wire  stress.  U  is  the  effective  mass  moment  of  inertia 
of  the  arm,  gripper,  and  the  payload,  and  c  is  the  torsional 
damping  coefficient  approximating  the  net  joint  friction.  For 
simplicity,  the  above  dynamic  equation  can  be  expressed  as 

Ie9  +  h{9,9)  =  Tw{cr) 


(2) 
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where  h{9, 9)  now  includes  viscous  damping  and  the  spring 
and  gravitational  terms.  The  SMA  wire  strain  rate  e  and  joint 
angular  velocity  9  are  related  kinematically  as 


where  Tp  is  pulleys  radius  and  Iq  is  the  initial  length  of  SMA 
wire. 


Fig.  2.  Block  diagram  of  control  design  including  both  feedforward  and 
feedback  loops 


B.  Thermomechanical  Model 

As  detailed  in  the  work  of  Liang  and  Rogers  [14],  the 
thermomechanical  behavior  of  SMAs  can  be  described  in 
terms  of  strain  (e),  martensite  fraction  (0,  and  temperature  (T). 
In  the  most  general  form,  the  mechanical  constitutive  equation 
is 

da  =  D{e,  0  T)de  +  f](e,  0  T)d^  +  0(e,  0  T)dT  (4) 

where  D{e,^,T)  is  representative  of  the  modulus  of  the 
SMA  material,  fl(e,0r)  is  the  transformational  tensor,  and 
0(e,0T)  is  related  to  the  thermal  coefficient  of  expansion. 

C.  Heat  Transfer  Model 

The  assumed  SMA  wire  heat  transfer  equation  consists  of 
electrical  heating  and  natural  convection: 

dT 

mcp— =  —  -  hAc{T  -  Too)  (5) 

where  R  is  resistance  per  unit  length,  Cp  is  the  specific  heat, 
m  is  mass  per  unit  length  and  Ac  is  circumferential  area  of  the 
SMA  wire.  Also,  V  is  the  applied  voltage.  Too  is  the  ambient 
temperature,  and  h  is  heat  convection  coefficient. 

III.  Control  Design 

In  order  to  create  an  improved  controller  design  we  first 
attempt  to  produce  a  linear  control  law  based  on  stress  rather 
than  angular  displacement.  A  simple  proportional-integral  (PI) 
controller  is  then  used  to  input  a  voltage  to  the  wire  by 
tracking  the  desired  stress.  It  is  worth  noting  that  much  of  the 
discussion  provided  here  is  simplified  somewhat  for  breadth. 
A  more  thorough  discussion  of  variable  structure  control  of 
robotic  manipulators  can  be  found  in  [15]-[18]. 

Recalling  the  dynamic  equation  of  motion  for  the  SMA 
manipulator.  Equation  2,  we  choose  the  input  to  be 

T^(a)  =  IcV  +  h(9,9)  (6) 

which  cancels  the  effects  of  damping  and  the  torque  due  to 
both  gravity  and  the  spring.  This  results  in  the  linear  equation 

9  =  v  (7) 

The  task  is  then  to  design  v  in  a  way  that  the  tracking  error 
converges  to  zero.  Defining  the  error,  e,  as  the  difference 
between  the  desired  position,  9'^,  and  the  true  position,  we 
can  define 

V  =  9‘^  —  kie  —  koe  +  Ksat{s)  (8) 


where  kg  and  ki  are  the  position  and  velocity  error  gains  and 
the  saturation  term  is  added  to  account  for  model  uncertainty. 
Combining  Equations  7  and  8  results  in 

e  +  kiC  +  kgc  =  K  sat{s)  (9) 

which,  for  positive  fco  and  ki,  guarantees  that  the  error  will 
converge  to  zero  in  finite  time  [17].  The  error  dynamics  can 
thus  be  controlled  by  a  process  of  selecting  the  the  error  gains, 
ko  and  ki,  and  the  sliding  control  gain,  K. 

Now,  given  a  desired  arm  angle,  the  desired  torque  can  be 
calculated  by  Equation  6.  Note  that  this  desired  torque  can  be 
converted  to  a  desired  stress  via  the  simple  kinematic  relation 


Alsin9yj 

where  A  is  the  cross-sectional  area  of  the  wire,  I  is  the  length 
of  the  arm,  and  9cu  is  the  constant  angle  between  the  SMA 
wire  and  the  arm.  The  problem  of  controlling  the  SMA  arm, 
however,  still  remains;  that  is,  how  to  derive  an  appropriate 
voltage  given  the  desired  torque.  This  problem  is  solved  using 
a  proportional-integral  (PI)  controller.  As  shown  in  the  block 
diagram  of  Eigure  2,  given  the  actual  and  desired  stress,  the 
variable  structure  controller  calculates  the  appropriate  stress. 
It  is  assumed  that  the  arm  angle  can  be  measured  with 
an  encoder.  The  stress  in  the  wire,  however,  can  not  be 
readily  measured.  To  alleviate  this  problem  a  Kalman  filter  is 
employed  to  find  the  stress  based  on  angular  position  and  the 
voltage  input.  The  design  of  the  Kalman  filter  is  not  discussed 
here,  however  the  reader  is  directed  to  [19],  [20]  for  a  detailed 
discussion.  Again  referring  to  Eigure  2,  the  predicted  stress 
from  the  Kalman  filter  is  compared  with  the  desired  stress  from 
the  variable  structure  controller  and  the  error  is  fed  through  a 
PI  controller  that  calculates  the  input  voltage  by 

V {f)  =  Kp{a  -  a'^)  +  Ki  f  (a  -  a‘^)  dt  (11) 

Jo 

where  Kp  and  Ki  are  the  proportional  and  integral  gain, 
respectively. 

IV.  Results 

The  following  sections  investigate  the  performance  of  the 
controller  design  presented  in  the  previous  section.  Numerical 
simulations,  performed  in  MATLAB/Simulink,  are  used  to 
implement  the  two-part  controller  design  with  the  purpose  of 
evaluating  tracking  performance  and  robustness.  A  linear  PID 
controller,  based  solely  on  angular  position  error,  is  used  for 
comparison. 
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Time  (s) 

Fig.  3.  Controller  performance  in  tracking  a  desired  step 


Fig.  4.  Tracking  desired  stress  for  a  45-degree  desired  angular  position 

A.  Control  Performance 

The  control  performance  in  the  tracking  of  a  step  input  can 
be  seen  in  Figure  3,  for  a  range  of  inputs  from  -10  to  90 
degrees.  Settling  times  of  around  2  seconds  are  shown  for 
lower  angles  and  around  3.5  seconds  for  larger  angles.  Also, 
some  overshoot  is  evident  for  angles  above  30  degrees.  This 
phenomena  is  due  to  the  sharp  reduction  in  gravity  induced 
stress  as  the  arm  passes  the  angle  of  maximum  torque  as  shown 
in  Figure  4. 

Figures  5  and  6  show  that  the  variable  structure  controller 
(VSC)  performs  better  than  a  PID  controller  based  solely  on 
the  angular  position  error,  e{t).  In  Figure  5  the  PID  gains 
were  optimized  at  22  degrees-representing  half  of  the  arm’s 
range  and  the  angle  of  maximum  stress.  It  is  evident  that  at 
angles  above  30  degrees,  PID  control  results  in  large  overshoot 
and  settling  time.  In  Figure  6,  the  PID  gains  were  optimized 
at  each  desired  arm  angle  to  demonstrate  that,  even  with  a 
variable  gain  PID  controller,  the  problems  with  overshoot  and 
settling  time  remain. 

In  Figure  7,  a  desired  angular  displacement  in  the  form 
of  a  40-degree  amplitude,  1  rad/sec  sinusoid  is  input  into 
the  controller.  After  approximately  2  seconds  relatively  close 
tracking  is  achieved.  The  only  significant  errors  occur  at  times 


Time  (sec) 

Fig.  5.  Control  performance  comparison  of  VSC  controller  and  a  simple 
PID  controller 


Fig.  6.  Control  performance  comparison  of  VSC  controller  and  a  simple 
PID  controller  with  optimized  gains  at  each  step 

when  the  angular  velocity  exceeds  the  desired  value.  It  is 
important  to  note  that  the  tracking  performance  of  a  sinusoidal 
signal-as  well  as  any  time-varying  input-is  dependent  on 
the  input  frequency.  There  are  time  constants  associated  with 
heating  and  cooling  of  the  wire  that  prohibit  high  frequency 
response. 

The  final  variable  input  tracking  simulation  is  based  on  a 
pseudo-random  user  input  of  Figure  8.  This  input  might  be 
typical  of  that  produced  by  a  joystick  control.  Again,  tracking 
performance  is  based  on  physical  limitations  of  the  SMA’s 
speed  of  response,  however  the  desired  path  is  followed  with 
good  accuracy. 

B.  Parameter  Uncertainty 

The  final  item  of  interest  is  robustness.  That  is,  given  some 
uncertainty  due  to  parameter  estimation  and  model  simpli¬ 
fication,  will  the  controller  continue  to  perform  adequately. 
In  the  following  simulations,  it  is  assumed  that  the  model 
is  accurately  defined,  however  we  include  some  parameter 
error  in  the  feedback  linearization  to  examine  the  effects  of 
uncertainty.  Specifically,  the  mass  moment  of  inertia,  /g,  is 
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Time  (s) 

Fig.  7.  VSC  controller  performance  in  tracking  a  40  degree  amplitude,  1 
rad/sec  sinusoidal  path 


Time  (s) 

Fig.  8.  VSC  controller  performance  in  tracking  random  user  commands 

varied  from  -50%  underestimation  to  h-50%  overestimation. 
Similarly,  damping  uncertainty  is  simulated  in  a  range  of  - 
50%  underestimation  to  +20%  overestimation.  The  simulations 
were  performed  by  inputting  the  false  estimations  into  the  VSC 
controller  of  Figure  2.  In  theory,  this  would  produce  error 
in  the  desired  stress  calculation  resulting  in  reduced  tracking 
performance. 

The  results  of  these  simulations  are  shown  in  Figures  9  and 
10.  Referring  to  Figure  9,  the  uncertainty  in  Ig  has  very  little 
affect  on  the  tracking  performance.  This  is  primarily  because 
the  inertia  term  is  tied  directly  to  sliding  surface  component  as 
described  by  Equation  6.  Conversely,  as  shown  in  Figure  10, 
both  overestimation  and  underestimation  of  the  true  damping 
results  in  either  large  overshoot  or  slow  settling  times.  This 
result  is  certainly  expected,  since  overestimation  in  damping 
would  result  in  a  lower  desired  stress,  and  underestimation 
would  result  in  a  higher  desired  stress. 

V.  Conclusion 

A  two-stage  stress  based  controller  was  designed  to  track 
desired  angular  positions  of  an  SMA  actuated  manipulator.  The 
controller  was  designed  such  that,  based  on  position  error,  a 


Fig.  9.  Effects  of  estimation  errors  in  the  mass  moment  of  inertia,  le',  errors 
are  varied  from  -50%  (underestimation)  to  -1-50%  (overestimation) 


Fig.  10.  Effects  of  estimation  errors  in  damping;  errors  are  varied  from  -50% 
(underestimation)  to  -1-20%  (overestimation) 


desired  wire  induced  stress  is  computed  by  a  variable  structure 
controller.  The  actual  stress  of  the  wire,  estimated  with  a 
Kalman  filter,  is  then  used  to  create  a  stress  error  signal. 
Based  on  this  error  signal,  a  proportional-integral  controller 
is  implemented  to  input  a  suitable  voltage  into  the  wire. 
The  benefit  of  the  additional  stress  feedback  loop  is  that 
the  nonlinearities  in  wire  stress  due  to  angular  position  can, 
in  effect,  be  cancelled.  Therefore,  the  problem  is  essentially 
reduced  to  a  hanging  mass  system. 

A  model  consisting  of  the  dynamic/kinematic  relations, 
the  SMA  constitutive  behavior,  and  heat  transfer  within  the 
wire  was  constructed  in  MATLAB/Simulink  for  the  pur¬ 
pose  of  evaluating  the  controller  design.  From  the  numerical 
simulations,  satisfactory  results  have  been  shown  in  terms 
of  tracking  both  stationary  and  time- varying  input  signals. 
Tracking  performance  of  time-varying  inputs,  however,  is 
limited  to  lower  frequency  signals  due  to  physical  limitations 
of  the  SMA’s  speed  of  response.  Furthermore,  the  stress-based 
tracking  method  is  shown  to  demonstrate  superior  performance 
over  a  simple  PID  controller-based  solely  on  angular  position- 
in  terms  of  rise  time,  settling  time,  and  overshoot.  With 
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regards  to  parameter  uncertainty,  it  was  shown  that  large 
errors  in  damping  estimation  can  result  in  undesirable  tracking 
performance. 
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